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Jacobi elliptic functions. We find that several of our Landen transformations have a rather different 
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1 Introduction 



Jacobi elliptic functions dn{x, m), cn(x, m) and sn(x, m), with elliptic modulus parameter m = k'^ [0 < 
m < 1) play an important role in describing periodic solutions of many linear and nonlinear differential 
equations of interest in diverse branches of engineering, physics and mathematics pp. The Jacobi 
elliptic functions are often defined with the help of the elliptic integral 

/ (1) 

Over two centuries ago, John Landen |2j studied the consquences of making a change to a different 



integration variable 



t 



{1 + k') z 



Vl - k^z^ 

This transformation yields another elliptic integral 

dt 



, k' = Vl - k^ = ^/T 



m 



[i + k')V{i-f'){i-i''f') 



1-k' 

1 + k' 



(2) 



(3) 



It readily follows that |3] 



dn 



cn 



il + k')u, 



{l + k')u, 



1-k 
1 + k' 



1-k' 
1 + k' 



sn 



(1 + k')u, 



1-k' 
1 + k' 



1 — (1 — k') sn^(ti,m) 
dn{u,m) ' 

_ 1 - (1 + A;') sn^(u, m) 
dn{u, m) ' 

(1 + k') sn(n, m) cn(n, m) 



(4) 

(5) 
(6) 



dn(ti, m) 

These celebrated relations are known as the quadratic Landen transformation formulas, or more simply, 
Landen transformations. They have the special property of providing a non-trivial connection between 
Jacobi elliptic functions involving two different unequal elliptic modulus parameters m and m, where 
the transformed modulus parameter is m = may be noted that in 

Eqs. Q) to ©, for < m < 1, m is always less than m and also lies in the range < m < 1. 

There are also similar formulas, known as the quadratic Gauss transformations where the 

transformed parameter m is always greater than m. These are obtained by a different change of 
integration variables 



t 



{l + k)z 
1 + kz^ 



(7) 
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This transformation yields the elhptic integral 

dt 



I 



I = 



2Vk 
1 + k 



(8) 



from which it follows that 



dn 



cn 



sn 



{l + k)u 

1 + k)u, 
{l + k)u 



Ak 



Ak 



Ak 
' (1 + A;)2 



1 — k sn^{u, m) 
1 + k sn2(ii, m) ' 

cn(n, m) dn(ti, m) 
1 + k sn2(n, m) ' 

(1 + k) sn(u, m) 
1 + k sn2(u, m) 



(9) 

(10) 
(11) 



In these Gauss transformation formulas, note that for any choice of m in the range < m < 1, the 
transformed modulus parameter m = 4-y/m/(l + -y/m)^ is always greater than m but lies in the range 
< 771 < 1. 

The above described Landen and Gauss transformation formulas are of order two. Subsequently, 
generalizations of these formulas to arbitrary order p have been studied [HI • The purpose of this paper 
is to obtain these generalized transformations by an entirely different method. In particular, we have 
recently shown [71|S| that a kind of superposition principle works for several nonlinear problems like Xcf)^ 
theory and for numerous nonlinear equations of physical interest [Korteweg deVries equation, modified 
Korteweg deVries equation, nonlinear Schrodinger equation, sine-Gordon equation]. Using the idea 
of superposition, we have obtained seemingly new periodic solutions of these nonlinear problems in 
terms of Jacobi elliptic functions. The purpose of this paper is to argue on mathematical as well 
as physical grounds that these solutions cannot really be new, but must be re-expressions of known 
periodic solutions j^. In the process of proving this, we discover the generalized Landen and Gauss 
transformations . 

As an illustration, let us focus on Eq. Q) first. Using the identity dn^(u, m) = 1 — m sn^(u, m), 
and changing variables to a; = (1 + k')u, one can re-write the order two "dn" Landen formula in 
the alternative form 



dn 



1-k' 
1 + k' 



a {di\[ax,m] + di\[ax + K{rn)^m]} , a 



1 



(1 + k') 



(12) 



Here, the right hand side contains the sum of two terms with arguments separated by K{m) = 
Jq^'^ d6[l — msin^ the complete elliptic integral of the first kind 01 [SJ. Our generalized Landen 
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formulas will have not two but p terms on the right hand side. We will show that the generalization 
of Eq. 0) [or equivalently Eq. (|12|) ]. valid for any integer p is given by 

p 

dn{x, rh) = a ^ dn[ax + 2{j — l)K{m) /p, m] , (13) 
i=i 

where ^ 

« = I E dn[2(i - l)K{m)/p, m] | , (14) 

and 

p 

m = (m - 2)a2 + 2a^^dn^[2(j - l)/C(m)/p,m] . (15) 
i=i 

However, if one starts from the "cn" or "sn" Landen formulas of order 2 as given by Eqs. © and 
the generalization to arbitrary order p is different depending on whether p is even or odd. The 
main results are: 

p 

cn{x,m) oc ^cn(5;j,m) , p odd, 

i=i 
p 

oc ^(— l)-'~^dn(xj, m) , p even, (16) 
i=i 

p 

sn{x, fh) oc ^sn(xj,m) , p odd, 

i=i 
p 

oc ]^sn(xj,m) , p even, (17) 

i=i 

where m is as given in Eq. (|15() . and we have used the notation 

Xj = ax + 2(j — l)K(m)/p , Xj = ax + 4(j — l)K{m)/p . (18) 

Note that all the above formulas have the same non-trivial scaling factor a of the argument x, as well as 
real shifts which are fractions of the periods of Jacobi elliptic functions. The richness of the generalized 
results is noteworthy and reflects the many different forms of periodic solutions for nonlinear equations 
which we have recently obtained [7| |H]. Some formulas involve the sum of p terms, the even p "cn" 
formula involves alternating + and — signs, and the even p "sn" formula has a product of p terms. In 
fact, there are also several interesting additional alternative forms for the above results which follow 
from use of several identities involving Jacobi elliptic functions which we have recently discovered 



4 



jlUlllllIT^ . For instance, we will later show that for odd p, the Landen formulas for dn, cn, sn which 
have been written above as the sum of p terms, can also be written as the product of the same p terms! 

Analogous to the generalized Landen transformations mentioned above, we also obtain generalized 
Gauss transformations in which the parameter m is greater than m. In particular, it is shown that 
Gauss transformations can be systematically obtained from the Landen transformations of arbitrary 
order p by considering shifts in the arguments by pure imaginary amounts and making use of the fact 
that the Jacobi elliptic functions are doubly periodic. 

By starting from the generalized Landen transformations discussed above, we show how to obtain 
Landen transformation formulas for other combinations of Jacobi elliptic functions like sn(x, m)cn(a:, m), 
sn{x,m)dn{x,m), cn(a;, m)dn(x, m), dn^(x,m), cn(x, m)dn(x, m)sn(x, m), etc. Besides, by integrating 
the Landen transformation formula for dn^(x, rh) we also obtain Landen transformations for the Jacobi 
zeta function Z(x,m) as well as for the complete elliptic integral of the second kind E[rh). Finally, by 
combining the recently discovered cyclic identities jlUlllllIT^ with the Landen transformation results, 
we obtain new types of Landen transformations. For example, a simple cyclic identity is jl2j 
p p 

^ dn^(xj, m)[dn(xj+r., m) + dn(xj_r, rn,)] = 2[ds(a, m) — ns(a, m)] ^ dn(xj, m) , (19) 
j=i i=i 

where a = 2rK{m)/p, Xj is as given by Eq. (|18|) and r is any integer less than p. Combining this 

identity with Eq. (|13|) . we immediately obtain a novel Landen formula relating dn(x, m) with the left 

hand side of identity (fT^ . 

The plan of this paper is as follows. In Sec. 2 we obtain solutions of the sine-Gordon equation 

in two different ways and by requiring these solutions to be the same, thereby obtain the Landen 

transformations for the basic Jacobi elliptic functions dn, cn, sn for both odd and even p. In Sec. 3 

we show that these Landen transformations can be expressed in an alternative form. Further, we also 

show that for a given p, the relationship between rh and m is the same, regardless of the type of Landen 

formula that one is considering. Using these results, in Sec. 4 we obtain Landen transformations for 

other combinations of Jacobi elliptic functions like sncn, sndn, etc. both when p is an odd or an 

even integer. In Sec. 5 we obtain Gauss transformations of arbitrary order in which the parameter rh 

is always greater than m. This is done by starting from the Landen formulas and considering shifts 

by pure imaginary amounts instead of real shifts. In Sec. 6 we combine the recently obtained cyclic 

identities with the Landen formulas, thereby obtaining some unusual new Landen transformations. 

Finally, some concluding remarks are given in Sec. 7. 
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2 Generalized Landen Formulas 



Given the diversity of the generahzed Landen formulas for Jacobi elHptic functions, it is necessary to 
estabhsh them one at a time. 

2.1 "dn" Landen Formulas 

To get an idea of our general approach, let us first focus on the proof of Eq. (|13|) in detail. Consider 
the periodic solutions of the static sine-Gordon field theory in one space and one time dimension, that 
is, consider the periodic solutions of the second order differential equation 

(pxx = sm<j) . (20) 

Note that the time dependent solutions are easily obtained from the static ones by Lorentz boosting. 
One of the simplest periodic solutions of Eq. (|2U|) is given by 

sin[(/)(x)/2] = dn{x + xo,rh) , (21) 

where xo is an arbitrary constant and this is immediately verifiable by direct substitution. It was shown 
in refs. j7||H] that a kind of linear superposition principle works even for such nonlinear equations as a 
consequence of several highly nontrivial, new identities satisfied by Jacobi elliptic functions |lU[llll[T^ . 
In particular, one can show _8, that for any integer p, one has static periodic solutions of Eq. H2U() 
given by 

p 

sin[^(x) /2] = Q ^ dn[ax + 2{j - l)K{m) /p, m] , (22) 

3=1 

where a is as given by Eq. H14|) . 

The question one would like to address here is whether solution 1)22(1 is completely new, or if it can 
be re-expressed in terms of simpler solutions like ()21|) . but where m and rh need not be the same. To 
that end, consider Eq. (|20|) more generally. We note that on integrating once, we obtain 

= C -2cos(j) , (23) 

where C is a constant of integration. Integrating again, one gets 

/ =x + xo, (24) 

J vC — 2cos(/5 
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where xq is a second constant of integration, which we put equal to zero without loss of generality 
since it corresponds to a choice of the origin of coordinates. On substituting sin(0/2) = -0, equation 
H24|) takes the form 

X . (25) 



Now the important point to note is that if we perform the integral for different values of C then we 
will get all the solutions of Eq. H2U|) . Further, if two solutions have the same value of C, then they 
must necessarily be the same. As far as the integral (|25|) is concerned, it is easily checked that the 
three simplest solutions covering the entire allowed range of C are 

ip = sech X , C = 2 , (26) 

ip = dn{x, m) , C = 4m - 2 , (27) 

ip = cn{x/V^ , m) , C = 2 , (28) 

rh 

where < m < 1. Note that the constant C has been computed here by using Eq. H23() . which in 
terms of ip{x) takes the form 

C = 2 - 4^-' + 7^4^ • (29) 
1 — ip"^ 

Thus, whereas for the solution ()27() . C lies in the range —2 < C < 2, for the solution 1)281) . C lies 
between 2 and cxd. Note that for C < —2, there is no real solution to Eq. H25|) . 

Now the strategy is clear. We will take the solution (|22() and compute C for it and thereby try to 
relate it to one of the basic solutions as given by Eqs. (|26|) to ()28j) . One simple way of obtaining the 
constant C from Eq. ()29() is to evaluate it at a convenient value of x, say x = 0. In this way, we find 
that for the solution C is given by 

p 

C = -2 + 4a2(m-2) + 8a3^dn3(2(j-l)K(m)/p,m) . (30) 

Now, as m ^ 0, a = 1/p, dn(x,m = 0) = 1 and hence C = —2. On the other hand, as m — > 
1, K{m = 1) = oo, dn(x,m = 1) = sech x and hence a = 1 so that C = 2. Thus for solution (|22() . 
as m varies in the range < m < 1, the value of C varies in the range — 2 < C < 2. Hence it is 
clear that the solutions (|22|) and ()27|) must be same. On equating the two C values as given by Eqs. 
(|27|) and (|3U|) . we find that the two solutions are identical provided m and m are related by Eq. 1)15(1 
and hence the appropriate Landen transformation valid for any integer p is given by Eq. ()13j) . Note 
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that when p = 2, one recovers the Landen formula H12|). since dn(K{m),m) = k' , and m simplifies to 
(l-k'f/il + k'f. 



2.2 "cn" Landen Formulas 

Unlike the dn case, it turns out that in this case the Landen transformation formulas for odd and even 
p have very different forms. We first derive the form for the odd p case and then consider the even p 
case. 

As shown in ref. jS|, another periodic solution of the static sine-Gordon equation (|2Uj) is 



sin((/)(x)/2) = ai cn 



aix 



m 



+ 4(j — l)K{m) /p, m 



, p odd , 



with «! being given by 



-1 



ai = <^ ^cn[4(j - l)K{m)/p, 
Using Eq. (|29|) . we can now compute the corresponding value of the constant C. We obtain 



m 



(31) 



(32) 



(33) 



It is easily checked that since < m < 1, C varies from 2 to oo, and hence the solutions H28() and (|31() 
must be identical. On equating the two values of C as given by Eqs. (|28]) and (j33|l . we then find that 
for odd p, the Landen transformation is 



cn(3;, mi) = ai ^ cn 



m 



+ 4(j — l)ii'(m)/p, m 



(34) 



where mi is given by 



m 

mi = ^ 
af 



(1 - 2m) + 2mQi ^cn^[4(j - l)K{m)/p,m\ 



(35) 



What happens if p is an even integer? As shown in [SI , in that case another periodic solution of 
the static sine-Gordon Eq. is 

p 

sin(0(x)/2) = a2''^{—iy~^dii[a2X + 2{j — l)K{m)/p,m] , p even , (36) 



where 02 is given by 



«2 = \ ^(-1)^-Mn[2(j - l)K{m)/p,m] 



(37) 
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Using Eq. (|29() . the value of C for this solution is easily computed and we find that 2 < C < oo. On 
comparing with solution (|28|) we find that in this case the Landen formula is 

p 

cn{x, 1712) = a2'^i-iy~^dn[a2Vr^ X + 2{j - l)K{m)/p,rn\ , (38) 
i=i 

where 1112 is given by 

if ' y 

rh2 = — {im-2) + 2a2j2i-iy'^dn^[2{j - l)K{m)/p,m]\ . (39) 

«2 [ ,=1 J 

As expected, in the special case of p = 2, we immediately recover the Landen formula 
2.3 "sn" Landen Formulas 

As in the cn case, here too the Landen formulas for even and odd values of p have very different forms 
and so we consider them separately. 

We start from the sine-Gordon field equation 

4>xx - (ptt = sin (/) , (40) 

and look for time-dependent, traveling wave solutions with velocity v > 1 (which are called optical 
soliton solutions in the context of condensed matter physics). In terms of the variable 

Eq. H40|) takes the simpler form 

(pr,ri = - sin (/) . (42) 

Note that the only change from Eq. H2U|) is an additional negative sign on the right hand side. On 
integrating this equation once, we obtain 

</,2 = + 2 cos (/> . (43) 

On integrating further, we get 

V + Vo , (44) 



^/C~+Ycos'< 

where r/o is a constant of integration which we put equal to zero without loss of generality. Substituting 
sin(^/2) = ip, yields 



If we now perform the integral for different values of C, then we get all the solutions. It is easily 
checked that the three simplest solutions of Eq. (|45j) covering the entire allowed range of C are 

il) = tanhr/ , C = 2 , (46) 

-(^ = "v/jfi sn (ry, m) , C = 4m — 2 , (47) 
tp = sn{-!^,m) , C = — -2, (48) 



where < m < 1. Note that the constant C has been computed here using Eq. (|43|) . which in terms 
of '0 takes the form 

C = -2 + 4^2 + . (49) 

1 — 

Thus, for solution ()47() . C is in the range —2 < C < 2, whereas for solution H48|) . C lies between 2 and 
oo. Note that for C < —2, there is no real solution to Eq. (|45jl . 

Using appropriate linear superposition, it was shown in ref. [H] that for odd p one of the solutions 
of Eq. is given by 

p 

sin{(p{r])/2) = ^/ma''^sn[ar| + A{j — l)K{m) /p,m] , p odd , (50) 
i=i 

with a being given by Eq. H14|) . Using Eq. (|49() . we can now compute the corresponding value of C. 
We find 

C = -2+'-^, (51) 
af 

where a, ai are given by Eqs. H14|) and (|32|1 respectively. It is easily checked that since < m < 1, C 
has values between -2 and 2. Hence the solutions H47|) and H5U|) must be identical. On equating the 
two values of C as given by Eqs. ()47() and H51|). we then find that the "sn" Landen transformation 
formula for odd p is given by 

p 

sn(x, 7713) = ai ^ sn[ax + 4{j — l)K{m) /p, m] , (52) 
with 7713 and m being related by 

7713 = . (53) 

af 

Finally, we turn to the "sn" Landen transformation formula for the case when p is an even integer. 
One can show |H] that in this case, a solution to Eq. 1)42(1 is given by 

p 

sm{(j){r]) /2) = m^^'^aAQ Y[ sn[ar] + 2{j ~ X^Kijn) /p, tti] , p even , (54) 
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with a being given by Eq. H14|) and Aq defined by 

p-i 

= sn{2jK{m)/p,m) . (55) 
i=i 

Using Eq. (|49|) . we can now compute the corresponding value of C. We obtain 

C = -2 + AmPa^A^ . (56) 

It is easily checked that since < m < 1, the value of C varies between -2 and 2 and hence the 
solutions H47() and (|54() must be identical. On equating the two values of C as given by Eqs. ()47() and 
H56|). we find that for even p, the Landen transformation formula is 

p 

Aq a sn(x, 7714) = sn[aj; + 2{j — l)K{'m) /p, m] , (57) 
i=i 

with m4 given by 

m4 = rrfa^Al . (58) 

Not surprisingly, for p = 2 we recover the Landen transformation formula It is amusing to notice 
that as m — > 0, 

p-i 

Ao{p, m = 0) = Jl sin(i7r/p) = . (59) 
i=i 

At this point, we have generalized all three of the celebrated two hundred year old p = 2 Landen 
formulas [Eqs. @, ©, to arbitrary values of p, the generalization being different depending on 
whether p is an even or odd integer. In the next section, we re-cast the formulas in even simpler 
form. Although several of these Landen formulas are already known p|, to our knowledge, they have 
never been derived via the novel approach of this article which makes use of solutions of nonlinear 
field equations. 



3 Alternative Forms for Landen Transformations 

3.1 Relation between the transformed modulus parameters rhi and m 

So far, we have obtained several seemingly different relationships between the transformed modulus 
parameters rh,rni,rn2,rn3,rn4 and m. Let us recall that while relation for m is valid for all p, 
relations (jSlJ and ()52() for mi and are only valid for odd p and relations (|38|) and H57j) for m2 
and m4 are only valid for even p. One would like to know if for any given p, the different Landen 
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transformations give the same relationship between rhi and m or not. For p = 2, it is known that 
rh = rn2 = [see Eqs. ((U to ©]• Similarly, with a little algebraic manipulation, one can show that 
for p = 3 the corresponding relations are 



where q = dn{2K(m)/3, m). Note that while deriving this result, use has been made of the fact that 



cn(4if(m)/3, m) = and that (7 satisfies the identity + 2(7^ — 2(1 — m)g— (1 — m) = 0. Similarly, 



using the relations dn{K{m) /2,rn) = dn(3i^(m)/2, m) = (1 — m)^/^ = t and dn(/C(m),m) = t^, it is 
easily proved that for p = 4, the relations are fh = m2 = rn^ = (1 — t)^/ (1 + t)^ . 

As p increases, the algebra becomes messier and it is not easy to write the corresponding relations 
between rhi and m in a neat closed form. However, we can still establish the equivalence of all rhi for 
any given p. For this purpose we equate the periods of the left and right hand sides of the various 
Landen transformation relations [^^, (jBH), (ESI), ^2^, (ISZI)]- We get: 

K{rh) = ^ , (61) 
pa 

_ K{m)^ 
A(mi) = ' (62) 



Kim) , , 

K{rh2) = ^-jd= , (63) 



i^(-3) = ^ , (64) 
pa 

Kim) , , 

KirhA = . (65) 

pa 

^Prom Eqs. H61() . ()64() and (|65() above, one immediately sees that fh = rh^ for odd p and rh = rh^ for 
even p. In order to establish that rhi = "^3 for odd p, we suitably rescale, square and add the Landen 
transformations (jHUl and (|52j) . This yields 



sn^(x, m) + cn^(-^^ — mi) = C , (66) 



where the constant C on the right hand side comes from making use of the cyclic identities |in| 111! I12j 
given in eqs. (|1U5() and 1)1061) . Since Eq. H66|) is valid for all and in particular for j; = 0, it follows 
that C = 1, and hence one must have 



sn (x,mj = sn ( ,mi). (67) 

ai\/rhi 
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This implies that "^'^5- = 1 and hence rfi-i = rh for odd p. Similar reasoning yields — %^ = 1 and 
m2 = ""T. for even p. 

Consequently, for odd p, the three Landen transformations are 

p 

dn(x, ?fi) = a ^ dn(xj, m) , (68) 
p 

cn(x,m) =ai^cn(xj,m) , (69) 
i=i 
p 

sn(x,m) = ai ^ sn(xj, m) , (70) 
i=i 

where a, ai are given by Eqs. (|14|) and (|H2|) . Xj,Xj are as given by Eq. (|18j) while m and m are related 
byEq. (fT5|) . 

Likewise, for even p, the three Landen formulas are 

p 

dn(x, TTi) = Q ^ dn(xj, m) , (71) 
p 

cn(x, m) = a2 ^(— l)-'^^dn(xj, m) , (72) 
i=i 

p 

^o"sn(x,m) = ]^sn(xj,m) , (73) 
i=i 

where a, 02,^0 are given by Eqs. (fTl|) . (pTTjl and respectively while fh and m are related by Eq. 

The numerical results for m as a function of m for various values of p ranging from 2 to 7 are 
shown in Table 1. Note that for any fixed value of p, as m increases from to 1, rh also increases 
monotonically from to 1 (but is always less than m). Also, for any given fixed value of m, m decreases 
monotonically as p increases and is always less than m. For this reason, Landen transformations are 
sometimes referred to as ascending Landen transformations [1]. 

3.2 Alternative Forms of Landen Transformation Formulas 

Recently we \U)\ I12j have obtained several new identities for Jacobi elliptic functions. Three of 
these, valid for odd p, are 

P (p-i)/2 ^ 

Y\_(in{xj,m) = Y[ cs^( , m) ^ dn(xj, m) , (74) 

j=l n=l P j=l 
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sn{x,m)= ^ ,„pli l^, r- (79) 



p (p-l)/2 p 

nsn(xj,m) = (-l/m)(P-^)/2 j^g2( , m) ^ sn(xj-, m) , (75) 

j=l n=l ^ j=l 

p (p-l)/2 p 

JJ^ cn(xj, m) = (1/m)*-*'^"'^^/^ ds^( , m) ^ cn(xj, m) . (76) 

j=l n=l ^ j=l 

Here ds(a;,m), etc. are defined by ds(x,m) = g^^^'^j ■ Hence, for odd the three Landen transfor- 
mation formulas can also be written as products (rather than sums) of p terms. In particular, for odd 
p the three Landen transformation formulas (|68() to H7U|) can also be written in the form 

1 ln=i dn(2nA [m) /p^ m) 
n?=i cn(xj, m) 

cn{x, rh) = ' -, (78) 

1 lra=i cn(4nA [m) /p, mj 

11^=1 sn(xj,m 

a(m) nn=i sn(4niir(m) /p, m) 

Similarly, for even p, on making use of the identity lllj 

(m)P/2 ]^sn(xj,m) = [ [| ns2(2ni^(m)/p, m)] ^(-l)^-^Z(xj, m), (80) 
i=i '^=l i=i 

where Z{u,m) is the Jacobi zeta function, we can rewrite the Landen formula ()73() for sn in the form 

p 

sn{x,m) = a2 ^(— l)-'^^Z(xj, m) . (81) 

We have not seen this particular form for the sn Landen transformation in the mathematics literature. 
Here the coefficient multiplying the right hand side of Eq. ()81() has been fixed by demanding consis- 
tency. From now onward, we shall mostly be using this form of the Landen formula rather than the 
one given by Eq. ((75|) . 

On comparing Eqs. (f75)) . (|5n|) and (|5T|) we obtain an interesting identity 

P-i f-i 



m 



p/2 



aa2 JJ^ sn(2jX(m)/p, m) = ns^{2nK{'m)/p,m) . (82) 

j=l n=l 

In concluding this section, we note that many of the above derived Landen transformation formulas 
lead to interesting known trigonometric relations by taking the limiting case m = m = 0. For instance, 
Eqs. dZSl), dZHl) and dZni) lead to 

p 

sinpx = 2^~^ Y[ sin[x + {j — l)vr/p] , p even, (83) 
i=i 
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cospx = 2^ ^ Y\_ cos[x + 2(j — l)vr/p] , p odd, 



sinpx = (— 4)~2~ II sin[a; + 2(j — l)'7r/p] , p odd. 



(84) 
(85) 



4 Landen Transformations for Products of Jacobi Elliptic Functions 

We shall now show that starting from the Landen formulas for the three basic Jacobi elliptic functions 
sn, cn, dn, we can obtain Landen formulas for their products and various other combinations. 



4.1 Any integer p 

We start from the basic Landen formula, Eq. (|68|) valid for any integer p. Differentiating it gives the 
Landen formula: 



sn(3;, m)cn(x, m) 



ma 



m 



2 P 



sn(xj, m)cn(a 



(86) 



For odd values of p, this can also be proved by multiplying the two Landen formulas [Eqs. (I^U]) and 
(fTO])] and using the cyclic identity 

p 



^ sn(xj, m)[cn(xj+r, m) + cn(xj_r, w-)] = , 



(87) 



obtained in refs. jlUl lllj . Here r = 1,2, ...,(p — l)/2. On the other hand, for even p, multiplying 
relations ()72|) and H81|) and comparing with identity ^86(1 yields the remarkable identity 



m ^ sn(xj , m)cn(xj , m) 



J2{-iy-'dn{x„m) 



Nontrivial results also follow from squaring any of the Landen formulas. For example, squaring 
the Landen formula and using the cyclic identity 
p 

^ dn(xj, m)dn(xj+r, ''7^) = p[dn(a, m) — cs(a, m)Z(a, m)] , a = 2rK{m)/p, 



yields the Landen formula for dn^, i.e. we get 



dn^(x, m) = 



Y,dn\xj,m) + 2Ad 



(90) 
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where 



Ad = ^ dn(xj, m)dn(xj, m) 
i<j=i 

{p-l)/2 

= p ^ [dn{2jK{m)/p,m) — cs{2jK{m)/p,rn)Z{2jK{m)/p)] , p odd , 

(p-2)/2 

= (p/2)Vl — m + p ^ [dn(2ji^'(m)/p, m) — cs(2ji^'(m)/p, m)Z(2jii'(m)/p)] , p even. (91) 

It is worth noting that the hnearly superposed solutions of KdV equation [Jj discovered recently, 
essentially correspond to this Landen transformation. 

Differentiating both sides of Eq. yields the Landen transformation for the product sn cn dn 
of the three basic Jacobi elliptic functions. We obtain 

v2 P 



dn(x, m)sn(x, m)cn(x, m) 



ma 



m 



dn(xj, m)sn(xj, m)cn(xj 



i, m . 



(92) 



On the other hand, if we integrate both sides of ()9U|) with respect to x, then we simulteneously 
obtain Landen transformations for both the Jacobi Zeta function as well as the complete elliptic 
integral of the second kind E{m). In particular, on using the well known formula ^[3] 



dn^(x, m) dx = Z(x, m) + 



E{m) 
Kim) 



(93) 



integration of Eq. ()9U|) yields 



Z(2;, rh) + 



E{m) 
Kim) 



-X = a X 



2Ad+p 



E{m) 
Kim) 



+ a j]Z( 



Xj,m) 



(94) 



Here the constant on the right hand side has been fixed by considering the equation at x = and 
using the fact that Z(x = 0,m) = 0, that it is an odd function of its argument x and that it is a 
periodic function with period 2K{m). It is now obvious that terms proportional to x on the left and 
right hand sides must cancel among themselves, since the remaining terms are oscillatory and do not 
grow. This immediately yields two remarkable Landen transformations: 

p 



Z(x, m) = a ^ Z{xj,m) , 



E{m) = a 



E{m) + 



2AdK{m) 



p 



(95) 
(96) 
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where use has been made of relation (|61|) . As expected, for p = 2 the Landen transformation agrees 
with the well known relation given in jlj. However, it seems that the generalized Landen transformation 
for Jacobi zeta functions is a new result. 

Landen formulas for higher order combinations like say dn"'{x,m) can be obtained recursively by 
differentiating the lower order Landen formulas ()86() and ()92p . For example, differentiating relation 
H86|) yields the Landen formula for dn^(2;,m) 

p p 



dn^(x, rh) = ^ dn^(xj, m) + a[2 — m — (2 — m)a^] ^ dn(a 



(97) 



4.2 Odd Integer Case 

We start from the two basic Landen formulas, Eqs. H69|) and ()7Up valid for any odd integer p. Differ- 
entiating and using the relation between rh and m, gives the following Landen formulas: 

p 



sn(x, m)dn(x, rh) = aai ^ sn{xj,m)dn{xj,m) 



p 



cn(x, m)dn(x, m) = aai cn{xj,m)dn{xj, 



(98) 



(99) 



It may be noted that these relations can also be proved by multiplying appropriate Landen formulas 
[Eqs. (|BS)) to ((70]) ] and using the cyclic identities analogous to (IH7|l obtained in refs. [inilTT] . 

Landen formulas for higher order combinations like sn^"^^(x,m) or cn^"+^(3;, m) can be obtained 
recursively from here by differentiating the lower order Landen formulas. For example, differentiation 
of relations (|98() and H99|) yield the following Landen formulas for sn^(a:,m) and cn^(x,m): 



cn^(x, rh) = af 



sn^(x, rh) = af 



cn"^ {xj , m) + 



■ p 



l-al 
a? 



sn^(xj, m) + 



+ 



1-Qf~ 



2a2 



^ sn(x'j, m) 



(100) 



(101) 



Before concluding the discussion about the Landen formulas for odd p, we want to point out the 
consistency conditions which we obtain by demanding dn^(x,m) + msn^(x,m) = 1, and sn^(x,m) + 
cn-^(x,m) = 1. In particular, using relations ()68|) to (jTU)) and demanding these constraints, yields 



1 m 

^ =p + 2{Ad + As) , = mp + 2{As + Ac) , 



at 



(102) 
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where in view of the cychc identities derived in |1U1 llll I12j we have 



(103) 



cn{2jK(m)/p) 



dn(2jK{m) /p, m)Z{2jK{m) /p, m) 



sn{2jK{m) jp) 



(104) 



V (p-l)/2 

= ^ cn(xj, m)cn(j;j, m) = p ^ 

i<j=l j=l 

and yld is as given by Eq. (|!JT|) . 
4.3 Even Integer Case 

We start from the two Landen formulas for even p given by Eqs. (|72() and H81|). Differentiating them 
and using the relation for m leads to the following Landen formulas: 



sn(x, m)dn(x, m) = maa2 ^( — 1)-' ^sn{xj,m)cn{xj, 



(105) 
(106) 



cn(x, m)dn(x, m) = 0:02 ^( — l)"* ^dn^(xj,i 

It may be noted that the relation (|106)) can also be derived by multiplying the Landen formulas 
((7T|) and (f7^ and using the cyclic identity 



dn(a;j, m)] ^(— 1)-' "'^dn(xj, m) = ^(— 1)-' "'^dn^(xj,m) 
j=i j=i j=i 



(107) 



On the other hand, multiplying relations (|71|) and H81j) and comparing with identity (jlOSj) yields a 
remarkable identity 



m^(— 1)-' "'^sn(xj, m)cn(xj, m) 



5:(-l)^-iZ(x„ 



^dn(j;j-,m) . 



(108) 



Landen formulas for higher powers like say sn^'^"^^(j;, m) can be obtained recursively from here by 
differentiating the lower order Landen formulas. For example, differentiation of the relations ()105() 
and ()106|) gives rise to the following Landen formulas for cn^(x,m) and sn'^(x,m): 



cn^(x, m) = 02 



^(-l)-'-Mn3(xj,m) + 



2a — ai m — 2 



m] 



(109) 



sn Ix, m ) = 02 



+ 02 
2a2 



^(-l)^-iZ(x,-, 



m)—ma\ y^( — 1)-^ ^sn(j;j, m)cn(xj, m)dn(xj, m) . (110) 
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Before finishing the discussion about the Landen formulas for even p, we want to point out the 
consistency conditions which we obtain by demanding dn^(a;,?TT-) + msn^(a;,m) = 1, and sn^(x,Tfi) + 
cn-^(x,m) = 1. In particular, using relations ()71|). 1)72(1 and ()81|) and demanding these constraints, 
yields 

^ -,2/ 



^ = dn^(xo,m) + 2Ad + 



^(-iy-iZ(x„ 



m] 



5:(-irMn(x, 



+ 



(111) 
(112) 



where is as given by Eq. (|9H1 . 



5 Gauss Transformation Formulas 

For p = 2, one has the Gauss quadratic transformation formulas as given by Eqs. @ to (|11|) in 
which m is greater than m. Gauss transformations are sometimes referred to as descending Landen 
transformations 4 . It is natural to ask whether Gauss transformation formulas can be generalized 
to arbitrary -p. In this context it may be noted that the Landen transformation generalizations which 
we have established so far are in terms of shifts involving the period K{m) on the real axis. We now 
show that the generalization of the Gauss transformation formulas to arbitrary order result from shifts 
involving the period iK'{m) on the imaginary axis. 

The procedure consists of starting with any Landen transformation formula, using the standard 
results 013 

dn(x,m') = dc{ix,m) , cn(x,m') = nc(ix,m) , sn(x,m') = —isc{ix,m) , (113) 

and then redefining ix = u. Note that m' = 1 — m while dc(x, m) = ^^f^, etc. In this way, we find 
that for odd p, the three Landen transformation formulas as given by Eqs. (|68|) to (|7U|1 take the form 

p 

dc(x, rh) = f3j2 dc[/3x + 2i{j - l)K'{m) /p, m] , (114) 
i=i 

p 

nc(x, m) = /?! ^ nc[/3x + 4i(j — l)i^'(m)/p, m] , (115) 
j=i 
p 

sc(x, m) = /?! ^ sc[/3x + 4i(j — l)K'{m) /p, m] , (116) 
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where ^ 

/3 = /?(m) = a(m' = 1 -m) = |^^dn[2(j - l)i^'(m)/p,m']| , (117) 

/3i = /3i(m) = ai(m' = 1 -m) = |^^cn[4(j - l)E:'(m)/p,m']| . (118) 

It might be noted here that the Landen formula (|114|) is in fact vahd for both even and odd p. On 
the other hand, for even p, instead of relations (|115|) and (|116j) we have the Landen formulas 

p 

nc(x,m) = l52^{-iy~^dc[l5x + 2i{j - l)K'{m)/p,m\ , (119) 
i=i 

p 

sc{px,m)(3Bo = '^^[^^ + ~ ^)K'{m)/p,m] , (120) 

where /32 and Bq are given by 

P2 = P2{m)=a2{m =l-m) = i^^{-iy-^dn[2{j-l)K\m)/p,m']^ , (121) 

p-i 

5o = Boim) = Ao{m') = J] sn[2jK'im)/p, m'] . (122) 

i=i 

For odd p, m is now given by 

>^ 

while for even p it is given by 

PI 

where /3, /9i,/?2 are as given by Eqs. (|117j) . ()118() and (|121() respectively. 

In the special case of p = 2, it is easily checked that rh = 4-y/m/(l + ^Jm)^ and as expected the 
relations (|114|) . (|119|) and (|12flj) reduce to the well known |^ Gauss transformation formulas as given 
by Eqs. © to (|11() . Note that there is an obvious, intersting relationship between the transformed 
modulus parameters mG{rri) and rai^{ra) in Gauss and Landen transformations respectively. For any 
choice of p, the Gauss transformation undoes the effects of the Landen transformation and vice versa. 
This means that mQ\mi^{m)\ = m and mi[?7iG(m)] = m. For p = 2, these relations are easily checked 
from Eqs. and Q. 

Before ending this section, it is worth remarking that just as we have considered Landen formulas 
where the shifts are in units of K{m) or iK'{m) on the real or imaginary axis respectively, we can also 
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l-m = (l-m)^, (123) 



1 - m = ^ , (124) 



consider Landen formulas corresponding to shifts in units of K{m) + iK'(m) in the complex plane. In 
this context it is worth noting that by starting from the Landen formulas for p = 2 as given by Eqs. 
Q to ©, changing m to 1/m, x to kx, and using the formulas 



dn(kx, — ) = cn(x, m) , cn(kx, — ) = dn(x, m) , sn(A;x, — ) = ksn(x, m) 
m m m 



one readily obtains the Landen transformations [31 



dn 



cn 



(k + ik')u, 



{k + ik')u, 



k-ik' 



k + ik' 



k - ik' 
k + ik' 



/X 2- 



sn 



(k + ik')u, 



k — ik 
k + ik' 



1 — k{k — ik') sn^(ti, m) 
cn(n, m) ' 

1 — k{k + ik') sn^(n, m) 
cn{u,m) ' 

{k + ik') sn{u,m) dn(ti,m) 



(125) 

(126) 

(127) 
(128) 



cn(u, m) 

The generalizations of these Landen transformations to arbitrary p is immediate. In particular, 
by changing ni{rh) to l/m(l/m), changing x to kx and using formulas (|125j) in the Landen formulas 
as given by Eqs. (|68j) to (|73j) . we obtain the corresponding Landen formulas for shifts by complex 
amounts in units of K(m) + iK'(m). For example, for any integer p, we get 

p 



cn(x, fh) = 6 ^ cn[5ia; + 2{j — l){K{m) + iK'{m)) /p, m] , 
while for odd p we have 

p 

dn(x, fh) = 62 ^ dn[52X + 4(j — l){K{m) + iK'{m)) /p, m] , 
i=i 

p 

sn(x, rh) = 61 ^ sn[62X + 4(j — l){K{m) + iK' {m))/p, m] , 
i=i 

with m and m being related by m = 61/6'^. 

On the other hand, the corresponding Landen formulas for even p are 

p 

dn(x,?fi) = 62 y^(— l)-^~^cn((52x/-v/m + 2{j — l){K{m) + iK' {m))/p,'m) , 
i=i 

p 

52-Dosn(x, rh) = ^Jm JJ^ ■sn.{b2X j \pm ^ 2(j — \)(K{m) + iK' {m))/p^ m) , 
i=i 

where m is given by fh = 62/6'^- Here, 6, 5i, 82 are given by 

p 

(5^cn[2(j - l){K{m) + iK'{m))/p,m] = 1 , 



(129) 

(130) 
(131) 



(132) 
(133) 

(134) 
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V 

(5i^dn[4(j - l){K{m)+iK'{m))/p,m] = 1 , (135) 
p 

'^2l](-lF"^cn[2(j - l){K{m) + iK' {m)) /p,m] = 1 , (136) 
i=i 

while -Do is given by 

P-i 

Do = J]^ SYi[2j{K{m) + ii^'(m))/^),m] . (137) 
i=i 

It is easily checked that for p = 2 the Landen formulas H129() , 1)132^ and H133() reduce to the well known 
"complex" Landen formulas (|127() . H126() and H128() respectively. 



6 Landen Transformation Formulas and Cyclic Identities 

Recently, we jlUl lllj have obtained a large number of cyclic identities where combinations of Jacobi 
elliptic functions at different points are expressed in terms of sums like J2^=i dn{x+2{j—l)K{m) /p, m), 
^^^-i^(— l)-'^"'^dn(a; + 2(j — \)K{m) /p,m), etc. Now the remarkable thing is that it is precisely these 
sums for which Landen and others have obtained the famous transformation formulas mentioned above. 
Thus by combining the cyclic identities and the Landen formulas given above, we can obtain a wide 
class of generalized Landen transformations for many combinations of Jacobi elliptic functions. 
For example, a simple cyclic identity valid for both even and odd p is 

p p 

dn^(xj, m)[dn(x.,-+r; ^) + dn(a:j_r, m)] = 2[ds(a, m)ns(a, m) — cs^(a, m)] ^ dn(j;j, m) , (138) 
i=i i=i 

where Xj is given by Eq. (|18|) and a = 2rK{m) /p. On combining with relation (|(i8|) . we obtain the 

Landen formula for the left hand side of (|138|) given by 

p 

dn^(xj, m)[dn(x.,+r, m) +dn(xj_r,m)] = (2/a)[ds(a, m)ns(a, m) — cs^(a, m)]dn(x, m) , (139) 

i=i 

where fh and m are related by Eq. ()15() . 

Proceeding in this way, using the many identities obtained by us |in[ 111! I12j , we obtain a huge 
class of new Landen transformation formulas. As illustrations, we give four examples: 

p 

m ^ cn(j;j, m)[sn(xj-|_r, "z) — sn(xj_r, w-)] = (2/a)[ns(a, m) — ds(a, m)]dn(x, rfi) ; (140) 
p 

^ dn^(xj, ?7i)[dn(xj-|-r7 ^) — dn(xj_,.,m)] = — (2m/A)cs(a, m)cn(j;, m)sn(x, rfi) , (141) 
i=i 
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where A is af or respectively depending on whether p is an odd or an even integer, and a, ai, 02 
are given by Eqs. (fTl|) . (|32|) and p7|) respectively; 



^ dn^(xj, m)dn^(xj+r., m) = — (2/a^)cs^(a, m)dn^(2;, m) 

+ 4^dcs^(a, m) + j5[cs^(a, m) + ds^(a, m) — 2cs(a, m)ds(a, m)ns(a, m)Z(a, m)] , (142) 



where A^i is given by Eq. H91() : 



dn^ (xj , m) [dn(x-,-|-r ; ^) — dn(xj_r.,m)] = — (2m/yl)cs(a, m)cn(a;, m)sn(x, m)dn(x, m) , (143) 

i=i 

where A is aaf for p odd and maa^ foi' P even. 

We now present four examples of novel Landen formulas which are only valid for odd p: 

p 

m''^sn'^{xj,rn)[sn{xjj^r,^) + sn{xj-r,m)] = — (2/ai)[ds(6, m)cs(6, m) — ns^(6, m)]sn(x, m) ; (144) 
i=i 

p 

rri' ^ sn^(xj, m)[sn^(xj+r, m) — sn^(xj_r, "i)] 

= (2/Qai)ns(6, m)[2ds(6, ?Ti)cs(6, m) + ns^(6, m)]cn(x, m)dn(2;, m) ; (145) 

p 

m ^ cn^(xj, m)[cn(x_,+r, ""T") + cn(xj_r, m)] = (2/Qi)[ns(6, m)cs(6, m) — ds^(6, m)]cn(x, m) ; (146) 
i=i 

p 

m ^ cn(xj, m)sn(xj, m)dn(xj, m)[sn(xj+r, m)dn(xj+r, m) — sn(xj„r; ?7i)dn(xj_r, m)] 

= (2/aai)ds(6, m)[ns(6, m)cs(6, m) + ns^(6, m) + cs^(6, m)]sn(x, m)dn(x, m) , (147) 

where Xj,Xj are given by Eq. ()18|) and b = 4rK{m)/p. 

On the other hand, here are five examples of novel Landen formulas which are only valid for even 

p and odd r < p: 

p 

m^(— l)-'~"^sn(xj,m)[cn(xj+r,m) — cn{xj-r,'m)] = (2/a2)[ns(a, m) + ds(a, m)]cn(x/a, m) ; (148) 
p 

m ^(— l)-'^^sn(xj, m)dn(xj, m)[cn(xj+r, 'm)dn(xj+r; i^) + cn(xj_r, m)dn(xj_r, m)] 
i=i 

= — (2/aa2)cs(a, m)[ds(a, m) — ns(a, m)]sn(x/a, m)dn(x/a, m) ; (149) 
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p 

^(— l)''~"'^dn(xj, m)dn(xj+r., m) = — (2/a2)cs(a, m)sn(x/a, ?fi) ; (150) 
i=i 

p 

^(— l)-'~^dn(xj, m)dn(xj_(_r; ?7i)dn(xj+2r) '"T')dn(xj+3r, m) 
i=i 

= (2/a2)[cs(a, m)cs(2a, m)cs(3a, m) + cs^(a, m)cs(2a, m)]sn(x/a, m) ; (151) 

p 

y^(— l)-^~^dn^(xj, m)[dn(x;/+r; m)+dn(x-,-r, m)] = (2/aa2)ns(a)ds(a)cn(x/a, m)dn(x/a, m) . (152) 
i=i 

7 Conclusion 

The results of this paper clarify the relationship between the well known periodic solutions of various 
nonlinear differential equations and those obtained recently by us using the idea of judicious linear 
superposition [71 |H]. In fact, we obtain a deep connection between the highly nonlinear Landen 
transformation formulas involving changes of the modulus parameter and certain linear superpositions 
of an arbitrary number of Jacobi elliptic functions. Further, we have found that using the identities 
with and without alternating signs along with Landen transformations, one can obtain novel Landen 
formulas for suitable combinations of Jacobi elliptic functions. In this context, it is worth recalling 
that recently we have also obtained cyclic identities with arbitrary weight factors and a large number 
of local identities T7. It would be very interesting if these could also be profitably combined with the 
Landen transformations. 
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Table 1: A table showing the transformed modulus parameter fh appearing in the generalized Landen 
transformation formulas as a function of the modulus parameter m and the number of terms p in the 
formula. The values of rh have been computed using Eq. H15|) . As discussed in the text, the same 
values of m could equally well have been obtained from the equivalent expressions Eqs. or H53() 
for odd integers p, and Eqs. H39|) or H58|) for even integers p. Note that the table also gives the change 
of modulus parameter for Gauss transformations if one interchanges the roles of m and rh. 



m 


rh{p = 2) 


rh(p = 3) 




m(p = 4) 




m{p = 5) 




rh{p = 6) 




rh{p = 7) 



































0.25 


.5155 X 10-2 


.9288 X 10" 


■4 


.1669 X 10" 


■5 


.3000 X 10- 


■7 


.5392 X 10- 


■9 


.9693 X 10" 


11 


0.5 


.2944 X 10"^ 


.1290 X 10" 


2 


.5580 X 10" 


■4 


.2411 X 10" 


5 


.1042 X 10" 


6 


.4503 X 10" 


-8 


0.75 


.1111 


.1005 X 10" 


-1 


.8666 X 10-- 


-3 


.7438 X 10- 


A 


.6381 X 10- 


-5 


.5475 X 10" 


-6 


0.9 


.2699 


.4311 X 10^^ 


-1 


.6158 X 10" 


-2 


.8655 X 10- 


3 


.1213 X 10" 


3 


.1701 X 10" 


-4 


0.99 


.6694 


.2506 




.7283 X 10- 


-1 


.1963 X 10" 


1 


.5185 X 10" 


-2 


.1362 X 10" 


-2 


0.999 


.8811 


.5292 




.2374 




.9312 X 10" 


1 


.3464 X 10" 


1 


.1264 X 10" 


-1 


0.9999 


.9608 


.7446 




.4481 




.2293 




.1080 




.4891 X 10" 


-1 


0.99999 


.9874 


.8721 




.6374 




.3973 




.2239 




.1193 




1 


1 


1 




1 




1 




1 




1 
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